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ABSTRACT 


A  rotating  ensemble  of  bodies  with  angular  periodicity  rotates  in  the 
presence  of  a  similar,  stationary  ensemble,  of  generally  different 
periodicity.  The  structure  is  illuminated  from  an  electromagnetic  source. 

The  scattered  field  is  modulated  periodically  as  the  scatterer  rotates,  and 
contains  a  discrete  spectrum  of  frequency  components. 

The  scattered  spectrum  is  predicted  through  electromagnetic  field  theory. 
The  theory  has  been  developed  such  as  to  exploit  the  angular  periodicities  of 
the  ensembles  and  thereby  reduce  the  computational  load  by  a  considerable 
factor. 

The  spectrum  consists  of  lines  separated  by  the  "body  rate",  i.e.  the  rate  of 
rotation  multiplied  by  the  number  of  bodies  in  the  rotating  ensemble.  The 
total  bandwidth  is  several  times  greater  than  that  for  the  rotating  ensemble 
in  free-space,  due  to  electromagnetic  interaction  between  the  two  ensembles. 


1. 


Introduction 


The  radar  scatter  from  jet  engines  is  periodically  modulated  by  their  rotation, 
resulting  in  a  line  spectrum  with  features  characteristic  of  the  engine.  Two 
previous  studies  [1,2]  of  the  radar  scatter  from  rotating  bodies  with  angular 
periodicity  have  been  undertaken  towards  predicting  the  spectrum. 

In  the  first  study,  a  single  stage  of  the  engine  compressor  was  modelled  as  an 
ensemble  of  coplanar  wires  radiating  from  an  axis.  Conclusions  were  drawn 
concerning  the  width  of  the  spectrum  and  the  line  separation,  and 
predictions  were  attempted  of  the  effects  of  other  features  of  engine 
compressors,  such  as  fan  blades  rather  than  wires,  multiple  stages,  the 
engine  intake,  stator  stages  etc. 

In  the  second  study,  a  limitation  of  the  wire  model  was  addressed;  the  wires 
were  replaced  by  an  ensemble  of  bodies  with  arbitrary  shape  but  angular 
periodicity.  This  extension  confirmed  that  many  of  the  earlier  conclusions 
remained  valid  for  a  single  stage  of  the  engine,  i.e.  for  a  set  of  blades 
mounted  on  a  central  hub. 

The  present  work  extends  the  model  further  by  including  a  stationary 
structure  in  the  presence  of  which  the  original  structure  rotates.  The 

stationary  structure  is  an  ensemble  of  conducting  bodies  with  angular 
periodicity,  coaxial  with  the  rotating  ensemble.  The  stationary  ensemble 
may  be  used  to  model  one  set  of  stator  blades  in  the  engine,  or  a  circular 
engine  cowling,  or  both. 

In  both  the  earlier  studies  the  incident  field  was  represented  as  a  spacial 
spectrum  of  harmonics  of  the  azimuthal  angle  (j)  measured  about  the  z  axis. 
Through  this  device  it  is  found  necessary,  for  each  harmonic,  to  determine 
the  current  distribution  on  one  body  only,  for  one  azimuthal  angle  of 

incidence  at  one  time;  all  other  cases  are  related  to  the  first  through  a  phase 
shift.  This  current  is  the  solution  of  an  electric  field  integral  equation  [3]  for 
scattering  from  a  single  body. 

The  spacial  harmonic  technique  is  again  employed,  leading  to  a  large 
reduction  in  computing  storage  and  time. 

The  approach  taken  here  involves  the  development  of  a  dyadic  Green's 
function  for  the  stationary  ensemble.  This  builds  on  the  work  in  [2]  but 
requires  determination  of  the  scattered  field  in  the  more  general  case  where 
the  source  and  observation  points  are  in  the  neighbourhood  of  the  scatterer, 
whereas  in  the  earlier  work  these  points  were  remote. 

The  field  scattered  by  the  rotating  ensemble  is  determined  with  the  help  of 
the  Green's  function  for  the  stationary  ensemble.  The  field  incident,  in  the 
presence  of  the  stationary  ensemble,  on  the  rotating  ensemble,  is  written  in 
terms  of  this  Green's  function.  Periodicities  in  the  geometry  are  exploited  to 
reduce  the  computation  of  the  currents  on  the  rotating  bodies.  These 

currents  radiate  the  scattered  field  in  the  presence  of  the  stationary 
ensemble;  the  field  is  calculated  with  the  Green's  function  for  the  stationary 
ensemble. 

As  has  been  a  feature  of  the  previous  studies,  certain  key  properties  of  the 

scattered  field  are  apparent  from  the  mathematical  form  of  the  equations 
without  need  for  their  numerical  solution,  and  attention  is  drawn  to  these. 

Notably,  the  scattered  field  spectrum  consists  of  lines  separated  by  the  "body 
rate",  ie.  the  product  of  the  rotation  rate  and  the  number  of  bodies  in  the 
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rotating  ensemble.  No  new  frequency  components  are  introduced  by  the 
presence  of  the  stationary  ensemble.  There  is  an  extension  of  the  width  of 
the  spectrum,  beyond  the  limits  seen  for  the  body  rotating  in  free  space, 

which  is  attributable  to  new  forms  of  electromagnetic  interaction  between 

stationary  and  moving  bodies. 

The  results  presented  here  for  two  angularly  periodic  ensembles  of  bodies 
with  relative  rotation  are  valid  for  the  limiting  case  where  one  or  both 

ensembles  contain  only  one  body;  thus  in  particular,  the  field  scattered  from 
a  periodic  ensemble  rotating  in  the  presence  of  an  arbitrary  stationary  body 
will  have  a  spectrum  with  lines  separated  by  the  body  rate,  and  of  bandwidth 
greater  than  that  for  the  ensemble  rotating  in  free  space. 
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2.  The  stationary  ensemble 

In  this  section  is  described  the  geometry  of  the  ensemble  of  conducting 
bodies,  together  with  the  coordinate  systems  and  some  special  parameters  to 
be  used  below. 

An  ensemble  S  of  I  conducting  bodies,  with  surfaces  Sj,  i=0,l . I-l,  is 

arranged  about  an  axis  (in  the  coordinate  system  to  be  used,  the  z  axis),  see 
figure  2.1.  The  surfaces  are  the  similar  and  oriented  such  that  Si  is  generated 
by  rotating  Sq  about  the  z  axis  through  the  azimuthal  angle 


The  bodies  may  be  detached,  as  shown,  or  attached  to  one  another.  In  the 
latter  case  the  individual  bodies  must  be  defined  by  arbitrary  cuts  through 
the  common  region,  creating  an  ensemble  of  bodies.  This  is  illustrated  in 
figure  2.2  which  is  suggestive  of  the  stator  blades  of  an  engine,  fixed  to  the 
cowling;  the  dotted  lines  are  the  cuts. 

Rectangular,  cylindrical  and  spherical  coordinates  will  be  used  as 
convenient;  the  system  in  use  will  usually  be  apparent  from  the 
conventional  symbols:  (x,y,z),  (p,(t),z),  (r,0,(t)). 

If  a  point  on  surface  Sq  has  position  vector  and  coordinates 


with  the  conventional  relations 

Po  =  (Xq  +  ■  '*’0  =  ^  =  sin-‘ 

Pq  Pq 

then  a  point  of  position  vector 
f.  =(p.,(b.,z.)  =  (r.,0.,(j).) 

with 

p.=Po.  Xi=Zo.  rj=rQ,  9.  =  9^ 

lies  on  Sj. 


(2.2) 


(2.3) 


Points  related  like  Fq  and  r.  will  be  called  "corresponding  points".  Such 

points  are  illustrated  in  figure  2.3.  Wherever  Fq  and  F.  or  their  similarly 

subscripted  coordinates  appear  in  the  same  equation,  they  are  to  be 
understood  as  corresponding  in  this  sense. 

The  relation  between  Fq  and  Fj  may  be  written  in  the  form 
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r.  =  R- 
1  1  0 


(2.4) 


where  Rj  is  a  dyadic  which  rotates  a  vector  about  the  z  axis  through  angle 
O. .  In  rectangular  coordinates  (2.4)  corresponds  to  the  matrix  equation 


r. 

cosO. 

-sinO. 

0" 

Tn 

IX 

1 

1 

Ox 

r. 

— 

sin(t>. 

cosdb. 

0 

ly 

1 

1 

Oy 

r. 

0 

0 

1 

r/^ 

.  iz_ 

L  OzJ 

(2.5) 


where  ri„  etc  are  the  rectangular  components  of  r.. 

lA  \ 


Rj  may  be  expanded  in  terms  of  spherical  unit  vectors  thus: 


R.  =  f.f^  +  0.0^  +  0.^^ 
1  1  0  1  0  0 


(2.6) 


where  the  unit  vectors  subscripted  i  and  0  are  evaluated  at  f^  and  Tq 
respectively. 

It  is  convenient  to  define  a  more  general  rotational  dyadic  R(A(t))  which 
rotates  a  vector  about  the  z  axis  through  azimuthal  angle  A(|).  The  dyadic  Rj 
is  the  special  case  where  Ac})  =  <!>.. 

Note  that  the  dyadic  R  ^(A({)),  which  reverses  the  effect  of  R(A(t)),  is  related 
to  it  by 

R“^(A(t))  =  R(-A(l))  .  (2.7) 

We  remark  that  any  vector  function  of  position  c(f),  whose  scalar 
components  expressed  in  spherical  coordinates  are  not  functions  of  (j),  takes 
values  at  two  positions  fj  and  r2  related  by 


c(r2)  =  R(A(l)).c(rj) 

if  Fj  and  r2  are  related  by  r2  =  R(A<j)).r^. 


(2.8) 
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Figure  2.1:  The  ensemble  of  bodies  with  surfaces  Sq,  disposed  with 

angular  periodicity  about  the  z  axis.  A  plane  wave  is  incident  from  the 
direction  of  spherical  coordinates  0\i|)*.  The  scattered  field  is  observed  from 
the  direction 
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Figure  2.3:  "Corresponding"  points  Fq  and  r.  on  bodies  Sq  and  Sj.  The  vectors 

c^(r^)  and  c^(r. )  are  rotated  relatively  to  one  another  about  the  z  axis  in  the 
n  U  n  i  ■' 

same  way  as  Fq  and  F^. 
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3.  The  free-space  dyadic  Green's  function 

3.1  A  direct  derivation  from  the  dipole  field 

The  electric  field  E  at  a  point  of  position  vector  f  of  an  ideal  dipole  of 
current  moment  p  at  a  point  of  position  vector  F  may  be  written  in  terms  of 

the  free-space  dyadic  Green's  function  G^(f,F),  thus: 

E(f)-^^(f,F).p  (3.1) 

Through  manipulation  of  the  formula  for  the  field  of  an  ideal  dipole  [e.g.  3]  , 
G^(r,F)  is  found  to  take  the  form 


G^(r,F)  =  G’!  +  G"RR  (3.2) 


where 

G'(R)  =  r - ^(-l-jk.R  +  k^R^)e  ^  ,  (3.3) 

jcoe  4;tR^  ^  ^ 

11  9  9  ~j^A^ 

G"(R)  =  -^ - ^-(3  +  j3k.R-k^R^)e  0  ,  (3.4) 

jcoe  47tR^  ^  ^ 


R  =  r-F  (3.5) 

jco^t 

and  k^  is  the  wave  number  for  e  ^  time  dependence.  The  time  dependence 

is  here  suppressed,  as  is  conventional.  It  will  be  reintroduced  when 
convenient. 


The  field  of  a  current  density  distribution  J(r)  is  then 

E(r)  =  111  G°(r,F).J(F)dv' 

V 


(3.6) 


in  which  the  volume  V  is  that  containing  the  current. 
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3,2  A  Green's  function  in  terms  of  spherical  harmonics 
A  formula  for  G^(r,F)  as  a  summation  of  spherical  harmonics  is  given  by  [4]: 

^  (2-5  )  (n-m)!  ^ 

47t _ I _ n  ^  n(n  +  n  (n  +  m)!„  „ 


^  n  =  lm  =  0 

fo  n(n  +  l)(n  +  m)! 

m(4)  M^l)' 

e  e 

+  R(‘‘)  N<')'  . 
e  e 

r  >  F 

mn  mn 

mn  mn 

0  0 

0  0 

m'-*)  m<»' 

e  e 

n(‘)'  , 

e  e 

r<r' 

mn  mn 

mn  mn 

o  0 


where 


5  = 

fl,  m  =  0 

m 

[0  ,  m  0 

_  m 

(r)  =  + - 

e 

sin  9 

II  n  sin 


^(0  (f)  =  £(!L±i)z(i)pniCOSj^^j,p 


ICqI  n  n  sin 


t'i'^  ■  m  cos  -  m  m  sin 

^.^^(pm  in(|)0  + - p™  m<t)  (()) 

o  o  sin  sin9  o  cos 


(3.10) 


the  primes  on  these  functions  specify  that  they  are  to  be  evaluated  at  F , 

"®=jn(V>'4‘‘^  =  ''n  ■  <3.11) 

(the  spherical  Bessel  function  and  Hankel  function  of  the  second  kind  [5]), 


i®  =  JLAta(i)) 

»  kptar'n'’ 


(3.12) 


pm  ^pm(  0) 

n  n  ^  ’ 


(3.13) 


(the  associated  Legendre  function  [5])  and 
pm  _  n 

n  ae 


(3.14) 
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In  the  above,  the  formulas  of  [4]  have  been  changed  for  consistency 
present  notation  and  the  time  dependence  (vice 

If  the  substitutions 


Sin  md)  = - ,  cos  md)  = - 

2j  ^2 

are  made  in  (3.9)  and  (3.10)  and  the  results  are  substituted  into  (3.7), 
finds,  after  considerable  algebraic  manipulation 


G°(f,f)  = 


jk 


0 


2n  + 1 


4k 


n  =  Im  =  -n 


mn  -mn  nui  -mn 
mn  -mn  mn  -mn 


r  >r 
r  <  r' 


where 


j^(i)  (f)  ^  2(i)(j  JILpm  0  _  pm  ^)ejm(t) 
mn'  n  '■'sinG  n  n 


n(0  (r)  =  (£i!L±I2z(0pm  j  +  z(i)(pm  9  +  jJILpm^))  jmtj) 
mn'  '  '  k^r  n  n  n  '  n  ■’sinG  n 

Upon  substituting  (3.17)  and  (3.18)  into  (3.16),  and  using  the  relation 

p-m  _  (n  -m)!pni 
^  (n  +  m)!  ^ 

one  finds  a  further  expression  for  the  dyadic  Green's  function: 

G®(r,F)  =  -^  i  I 

where 


and  for  r  >  r', 


g 


rr 


n^(n  +  l)^ 


h  i  '  pmpm, 
n-^n  n  n 


with  the 

(3.15) 

one 

(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 


.  n-'n  n  n 

*^0 

+  ■•  ._rT^pmpm. 

^r({)'  J  .  n-'n  sin0'  n  n 

KqF 

n(^l^  .pTOpm, 

°0r  .  ,  n-’n  n  n 

KQr 


gQo.  =h  j  ' - - - P'^P^'+h  j  •p'^p"^' 

00  n-’n  sin  0 sin  0'  n  n  n-'n  n  n 


g  -_j(h  j  •JELpmpm.^h  j  .'_m_pmpm.) 
*0(t)  n-’n  sin0  n  n  n-’n  sin0'  n  n 

g  -j£(!±i2h  j  .J!Lpmpm, 

<!>r  I,  r’  ^  ^  sin0  n  n 

KQr 

gAfi.  =j(h  j  ’-^P™P"^’-Hh  j  •-ELpmpm.^ 
^d>0  n-'n  cinff  n  n  n-'n  cinfi  n  n  ^ 


g^^'  =h  j  'pmpm.^jij  j  . — m - pnipin. 

^({)({)  n-’n  n  n  n-’n  sin  Osin  0'  n  ” 


(3.22) 


while  for  r  <  r'  the  elements  g^, ,  g^,  etc  are  found  through  the  exchange  of 
hn  and  j„. 


(3.20)  may  be  further  rearranged  into  the  form 


G^(r,r)=  i 


m  =  -oo 


where 


(3.23) 


G®(r,r)  =  -^  I  _ll±L(il_E2:g  (rf) 

""  4k  j^^|;;i^Qn(n  +  l)(n-nn)!^nin' 


(3.24) 
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3.3  Remote  source 


(3.24)  might  be  specialised  for  this  case  by  substituting 
it  occurs  in  the  elements  of  g  see  (3.22),  the 

approximation 


for  h^^^(k;p,r') 
n  '  0 

large  argument 


where 


;n  +  1  S 


-jkor' 


kor' 


(3.25) 


A  simpler  approach  will  be  used.  It  is  based  on  the  fact  that,  when  the  source 
of  radiation  is  distant,  the  field  incident  on  the  ensemble  of  bodies  is  a  plane 
wave. 


Consider  a  point  source  at  the  remote  point  of  position  vector  f  =  (r' ,0' ,(j)' )  of 
vector  current  moment  p.  We  are  interested  in  the  field  of  the  dipole  near 

the  origin  and  so  consider  the  two  polarisations  p  =  0'  or  .  (The  third,  f' , 
produces  no  field  near  the  origin.)  The  incident  electric  field  at  the  origin  is 


E'(r  =  0)  =  pEQ 


where 


E.  =-pjk.Z.^ - 

0  0  0  4T±r' 


(3.26) 


(3.27) 


is  the  scalar  field  value  at  the  origin. 


The  field  in  the  vicinity  of  the  origin  is  essentially  a  linearly  polarised 
plane  wave.  Such  a  wave,  with  the  vector  value  given  in  (3.26),  has  electric 
field  vector  at  the  point  r  =  (r,0,(j))  in  the  vicinity  of  the  origin, 

jk„r(sin0sin0'  cos((t)  -  $' )  +  cos0cos0' ) 

E(r)  =  pEQe  ^  (3.28) 

We  consider  the  two  polarisations 


p  =  0'  =  f(sin0cos0'  cos((l)  -  $' )  -  cos  Osin  0' ) 

+  0(cos  0  cos  0'  cos((t)  -  (1)' )  +  sin  0 sin  0' )  -  ^  cos  0'  sin((t)  -  (])' ) 


(3.29) 


and 


p  =  ^’  =  -f  sin  0  sin((})  -  f )  +  0  cos  0  sinftj)  -  ({)' )  +  ^  cos((t)  -(])'). 


(3.30) 


Through  application  of  the  formulas  in  Appendix  A,  (3.28)  to  (3.30)  may  be 
manipulated  into  the  form 
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E‘(r)=EQ 

m  =  -oo 


(3.31) 


where,  for  p  =  9' , 


c^(f,0' )  =  c^(r,0' )  =  j^(-r(jsin9cos0'  j^(kQrsin0sin0'  )  +  cos9sin0' J^) 

-0(jcos0cos0' j  -sin9sin9'J  ) 
m  m 

.K  m  jk^rcos9cos0' 

+  (j)cos9' - J  )e  ^ 

kQrsin0sin0' 


and  for  p  =  , 


(3.32) 


0COS0J  )-^jj  )e 


jkQrcos0cos9' 


(3.33) 


in  which  J  is  the  Bessel  function  of  order  m,  J  is  its  derivative  with 
m  m 

respect  to  the  argument,  and  the  argument  in  all  occurrences  is 

kQrsin0sin0' ,  as  shown  explicitly  in  the  first. 

A  free-space  dyadic  Green's  function  for  a  distant  source  is  extracted  from 
these  results  in  the  following  fashion: 

With  the  basic  property  of  the  Green's  function,  Cf.  (3.1), 


E\r)  =  G®(f,F).p  . 


(3.34) 


(3.27)  and  (3.31)  to  (3.33)  lead  to 


G^(tf)  =  y  G®(rf) 
m  =  -oo 


(3.35) 


where 


)r'  ^ - (c®'  (r,9'  )0'  +c^'  (f,9'  w ) 

m'  ^r  -» oo  •’0  0  4irr'  ^ 


(3.36) 


It  will  sometimes  be  convenient  to  express  this  in  the  form,  (note  (3.27)), 


(3.37a) 
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This  equation  may  be  expressed  alternatively 

E, 


I  aP'c^  (r.e'). 


G®  (r,F)  . 

'r  — >  oo  p  r,,  ccini 

^  a,p 


0'  -0' 
where  C  is  the  component  of  in  the  a 
a;m  ^  m 


in  the  following  notation: 

a  =  r,e,({);  p'  =  0',(|)'  (3.37b) 


direction. 


3.4  Remote  field  point 

This  is  the  reciprocal  of  the  previous  case.  It  is  readily  found  that 
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G^(r,r)  =  i  G®(f,r) 

/f— >oo  ^  m  '^r— >oo 

m  =  -oo 


(3.38) 


where 


=  -jkQr 

>r  ~  .0)  +  (?  .0)) 

^  3'  =  r',e',0' 

(X,  P 


(3.39) 
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3.5  Green's  function  harmonics  and  rotation 

Consider  the  free-space  dyadic  Green's  function  given  in  (3.23)  and  (3,24) 
and  note  that  any  field  calculated  from  it,  see  (3.1),  involves  vectors  resulting 
from  the  operation  Expanding  this  in  terms  of  its  components  gives 


=  '''Sre’Pe'  +  Pr’  +899'P0'  +  % P*' 

+  *S(|,t'Pr'  +8,t,e’Pe'  +g$4,'P4,')* 


(3.40) 


in  which  p^, ,  Pg,  ,  p^, 
location  F . 


are  respectively  the  r',  9',  (})'  components  of  p  at  its 


It  is  seen  from  (3.22)  that  the  scalar  components  of  the  vector  in  the  right 
side  of  (3.40)  are  not  functions  of  (().  From  the  remark  preceding  (2.7)  it 
follows  that  the  vector  transforms  between  points  related  by 

f2  =  R(A(j)).rj  as  in  (2.7),  i.e., 


'  <3.41) 

From  this,  noting  that  p  is  arbitrary,  we  have  the  dyadic  identity 

Smn<'2’f>  =  R('3't’>-Imn<'r''>  ’  <3-«) 

and  hence,  noting  (3.24), 

(I(A(t)).r,F )  =  I(A(1)).U°  (r,r )  .  (3.43) 

This  property  will  be  shared  by  any  dyadic  whose  scalar  components, 
expressed  in  spherical  coordinates,  are  independent  of  (|). 

Similarly,  because  the  elements  of  G^(r,F)  are  independent  of  ({)', 

G^(r,I(A(})).F )  =  G^  (r,F  ).I~^(A(t))  .  (3.44) 


where  R  ^(A(())  is  the  dyadic  which  reverses  the  rotation  R(A(()). 

Finally,  it  is  convenient  here  to  note  the  transform  property  of  the  free- 
space  Green's  function 

^^(R(A(1)).  r,R(A(t)).F  )  =  R(A(l)).5^(r,F ).  R“^(A(})). 


(3.45) 
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4.  The  incident  electric  field  and  the  induced  currents 

We  now  consider  the  field  of  the  point  dipole  to  be  incident  on  the  ensemble 
of  bodies  described  in  Section  2.  To  distinguish  it  as  such  it  is  superscripted  i. 

(3.1)  and  (3.23)  give 

OO 

E’=  Z  (4.1) 

m  =  -OO 

where 


E*  (r)  =  GO(f.r).pej"'<'l’-''’' 


(4.2) 


and  G^(r,F)  is  given  by  (3.24)  or  (3.36)  as  appropriate. 

(4.2)  and  (3.43)  show  that  Ej^(r)  varies  between  corresponding  points  r^  and 


r.  as 

1 


(r.)=R..E^  (r„)e  ^ 


(4.3) 


in  which  O.,  defined  in  (2.1),  is  the  azimuthal  angle  between  the  points. 


(4.1)  and  (4.2)  represent  the  field  as  the  sum  of  a  series  of  harmonics  of  the 
azimuthal  coordinate  (|).  A  given  harmonic,  the  mth,  say,  takes  vector  values 
which  differ  at  the  corresponding  points  fj  and  Tq  on  the  bodies  Sj  and  Sq  by 

rotation  through  angle  about  the  z  axis,  and  phase  shift  The  rotation 

of  the  field  vector  is  the  same  as  that  of  the  bodies  themselves,  whence,  the 
field  experienced  by  each  body  relative  to  its  own  orientation  is  the  same 
except  for  a  phase  shift. 


It  is  concluded  that  the  vector  current  distributions  induced  by  the  mth  field 
harmonic  are  the  same  on  all  bodies,  except  for  rotation  and  phase  shift, 
thus : 


J. 

im 


^  1  1  Om 


(Io)e 


jmOj 


(4.4) 


where  Jjjjj(rj)  is  the  surface  current  density  at  on  Sj  induced  by  the  mth  - 

field  harmonic  E^  . 

m 


As  the  bodies  are  perfectly  conducting,  the  currents  flow  on  their  surfaces. 
The  quantities  in  (4.4)  are  surface  current  densities,  dimensioned  A/m.  The 

— c 

conventional  symbol  for  surface  current  density  is  J  ,  but  in  this  work  the 
superscript  s  will  be  suppressed. 


(4.4)  reveals  a  major  benefit  of  the  harmonic  representation  of  the  field:  it 
will  be  necessary  to  determine  the  current  distributions  on  only  one  of  the 
bodies,  since  all  are  similar  except  for  orientation  and  phase  shift. 
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5.  The  electric  field  integral  equation 

Each  current  harmonic  induced  on  the  surfaces  radiates  a  part  of  the 

scattered  field,  E^(r),  say.  For  the  present  we  shall  refer  to  this  as  the  mth 

"harmonic"  of  the  scattered  field,  although  this  will  later  be  seen  to  be  a 

misnomer.  The  mth  harmonic  of  the  total  field  is  E„(f)  plus  the  mth 

m  ^ 

harmonic  of  the  incident  field,  E'  (T).  The  total  mth  harmonic  satisfies  the 

m 

boundary  condition  that  everywhere  on  the  surfaces  Sq,...,Sj.j  its  tangential 

component  is  zero.  This  is  expressed  in  an  integral  equation  whose  solution 
is  the  mth  surface  current  density  distribution. 


The  field  E^(r)  is  produced  by  the  currents  on  all  the  bodies,  “ 

It  is  given  in  terms  of  the  dyadic  Green's  function  (Cf.  (3.6))  by  the 

sum 


i  =  0f.'onS. 
1  1 


I-l 


fO, 


(5.1) 


Where  dsj'  is  an  element  of  area  at  r.'  on  Sj. 

With  (2.4)  and  (4.4)  this  is  rewritten  in  the  form 


(5.2) 


i  =  0 


Tq  ouSq 


The  boundary  condition  is 

[E^  (r.)  +  E^  (r.)L=0.  all  r.  on  S.,  j  =  0,1,. ..,1-1  (5.3) 

j  j  T  J  J 

where  the  notation  [.]j  specifies  the  vector  component  tangential  to  the 
surface. 

With  (4.2)  and  (5.2),  (5.3)  becomes 


I-l  jmtl).  =n  =  =  _  = 

(Xe  >  II  a% 

i  =  0  ?  ’--C  J 


Tq’  onS^ 


=  -[G°  (f.,F).p]  e"  ''■> 
'■  m  j 


jm((t).  -4)’) 


(5.4) 


all  r.  on  S  ,  all  j. 
J  J 


In  the  right  side  of  this  equation  ‘s  determined  from  (3.24)  or 


(3.36)  as  appropriate. 
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(5.4)  must  be  satisfied  for  all  j.  Explicitly,  for  j  =  0, 

I-l  jmO.  _ 

[  S  e  ■_  J|G«(f(,,R,V).R,Jo„(Fo')dSo'], 

1  =  0  r^jOnSg 

=0  _  “  <t>’ ) 

^  ^  .allr^onS, 


‘0  0 


With  (2.4)  and  (4.4),  this  may  be  written,  after  some  rearrangement 
I-l  jmO.  .  =  =n  =  = 

^0«"So  (5.6) 


=rA  J***\  Y ;  T  / 

=  -tG°(rj,F).p]je  ],,  on 


jm(<t).-(j>') 


in  which  it  has  been  noted  that  O.  +  <I).  =  (I>.  .  and  R-  =R.  ^  R.  . 

1  J  i+j  1  J  i+j 

Using,  see  (3.45), 

Rj.CO{fo,f;').Rj-*=50(fj.fj^y)  (5.7) 

and  replacing  the  summation  index  i  by  i+j,  we  may  recover  (5.4)  for  r.  on  Sj. 
Thus  satisfying  the  condition  on  Sq  satisfies  it  on  all. 

In  (5.5),  multiply  both  sides  by  ^  and  change  the  order  of  summation  and 
integration  to  get 

I  It 


ro'onSo 


=0  _  _ 

=  0,  alltgOnSQ 


where 


=y  jm(I).=M  =  = 

G^;(fn.fn'>=  S  e  ‘GO(rg.R..fQ').Rj 


J  (f  ' )  =  J  (f  ' 


(5.10) 
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The  azimuthal  position  of  the  source  ({)'  does  not  appear  in  (5.8)  and  hence  the 
new  unknown  independent  of  (()’.  This  is  a  further  benefit  of  the 

present  approach:  if  we  are  interested,  as  we  shall  be,  in  source  positions 
with  a  range  of  azimuthal  angles  ())',  all  with  the  same  polar  angle  0',  it  is 
unnecessary  to  solve  more  than  one  integral  equation,  (5.8).  The  currents  on 
all  bodies,  for  all  azimuthal  incidence  directions  are  then,  from  (4.4)  and 
(5.10), 


J. 

im 


(r.)=  R..J 

I  1  nn 


(ro)e 


jm(<I).  -  {})' ) 


(5.11) 


A  case  of  special  interest  is  that  in  which  the  source  point  is  distant.  Let 

-3'  , 

.  3'  =  0'.<1>’  be  the  current  distribution  when  the  source  has 
m  0  r  ->  oo’  ^ 

_  "  -B'  _ 

polarisation  p  =  3'p-  Then  with  (3.37a)  and  (5.8),  ^be 

solution  of 


Tq'  ohSq 


_g'  _  jmt}). 


(5.12) 


all  Fq  on  Sq 
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6.  The  scattered  field 

The  scattered  field  is  the  sum  of  the  fields  radiated  by  all  the  harmonic 
currents  on  all  the  bodies.  It  is  given  by  the  formula 

—  OO  I  —  1 

E®(f)=  s  X  JJ  G“(f,fj')Jj,(rpds.',  ( 

^  = -ooi  =  Or.' onS. 

1  1 

Substitution  using  (3.23)  and  (5.11)  and  slight  rearrangement  leads  to 


I  =  — oom  =  — oo  i  =  0 

II  ““So 


To'onSo 


Finally  we  note  from  the  definition  of  Oj,  see  (2.1),  that  the  summation  on  i  is 
zero  unless  £  +  in  is  an  integer  multiple  of  I,  in  which  case  it  equals  I. 
Therefore  set  ^  =  nl  —  m  and  sum  over  integer  values  of  n  to  get 


m  =  -oo 


fo'onSo 


(6.3 
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7.  The  Green's  function  for  the  stationary  ensemble 

7.1  The  general  case 


’““S  — 

In  (6.3),  E  (r)  is  the  scattered  field  at  F  when  the  point  dipole  p  at  f  radiates 
in  the  presence  of  the  ensemble  of  bodies.  This  may  be  expressed  in  terms  of 
a  dyadic  Green's  function  relating  the  scattered  field  to  the  dipole  thus  (Cf. 
(3.1)): 


-ES(r)  =  G^(r,F).p 
P 


(7.1) 


G^(r,F)  (note  the  lower-case  superscript)  may  be  expanded  in  terms  of  its 
component  dyads  in  spherical  coordinates  as  follows: 

G^(r,F)=  I  G^f.,(r,F)dP'  a  =  r,0,(|);  3’=  F,9',(j)’  (7.2) 

a,P' 


Let  the  dipole  be  oriented  in  the  p'  direction,  i.e.  p  =  P'p;  for  this  case  let 
J^(ro')  be  the  solution  of  (5.8),  E^^  be  the  scattered  field  calculated  from 

sB* 

(6.3)  and  EJ^  be  its  a-component.  It  is  seen  from  (7.1)  and  (7.2)  that 

G®  .  =-E^P'(r) 

ap  p  a  ^ 


With  this  result  and  (6.3)  G^(r,F)  may  be  written  as  a  double  summation  of 

terms  G^„(r,F): 
mn'  ' 

G^(f,F)=  i  i 


m  =  — oon  =  — oo 


m  =  -oon  =  -oo  cx,B'  ^  ’ 

a  =  r,0,<|);  |3'=  r’,0',(t) 


(7.3) 


where  G^o,  is  the  aB’  element  of  G^„(r,F),  and 

ap  ;mn  ^  mn  v  » 


mn 
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(7.4a) 

in  which,  in  anticipation  of  later  developments,  I  has  been  replaced  by 
and  Sq  by  Sq. 


If  the  source  point  is  distant,  (7.4a)  should  be  evaluated  with  the  current 
density  in  the  integrand  given  by  the  solution  of  (5.12). 

If  the  observation  point  is  distant,  (7.4a)  becomes,  note  (3.39), 


fo'onSg 


1  c“,  (f'.e)/'  (f  'je  ^"“'’Ods  ■ 

|3Q'  =  rQ',eQ',<.Q' 1^0  °  |3o’;nlS-m  0  0 


(7.4b) 

The  total  field  produced  by  the  dipole  p  in  the  presence  of  the  ensemble  is 
the  sum  of  the  scattered  and  incident  fields,  see  figure  7.1.  The  dyadic 
Green's  function  for  the  total  field  is  the  sum  of  that  developed  in  Section  3 
for  the  free-space  case  and  that  in  (7.3).  Using  for  the  former  the 
representation  in  (3.23)  and  for  the  latter  (7.3)  we  have 


G^(r,F)=  I  X  (r,F) 


m  =  -oo 


where  G^(r,F)  (note  the  upper-case  superscript)  is  a  dyadic  Green's 
function  relating  the  total  field  to  the  dipole  radiating  in  the  presence  of  the 
ensemble  S,  and 

^m0(f-'’>  =  *n0G°(f,F)+5^0(f,F)  (7.6) 

where  5jj  is  the  Kronecker  delta,  equal  to  one  if  i  =  j,  otherwise  zero. 

(7.5)  and  (7.6),  with  the  definition  of  their  terms  in  (3.23)  and  (7.4), 
constitute  a  formal  solution  to  the  problem  of  the  dyadic  Green's  function  for 
the  ensemble  of  bodies. 
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7.2  The  dyadic  Green's  functions  and  rotation 

In  (7.4)  the  functions  of  (t)  or  0'.  It  follows  that 

1  llllf 

(r,F)  transforms  as  in  (3.43)  and  (3.44);  specifically, 


G^^(R(A(})).r,F)  =  R(A(})).G^^(r,F), 

=  =  =r  =  ,  (7-7 

G^(r,R(A({)).F)  =  G^(r,F).R-^(A({)) 

It  has  been  noted  in  (3.43)  and  (3.44)  that  G^(r,F)  transforms  in  the  same 

m 

“S  - 

fashion.  It  follows  from  (7.6)  that  G^^  (r,F)  also  transforms  in  this  way: 

mn  ■' 


G^(R(A<|i).f,F)  =  R(A(t.).G^(?,F), 
5^(r,f(A(|».F)  =  5^(r,F  ).f-‘(A(t,) 

These  two  formulas  together  give 

(r,F  ).R"l(A(t)) 


(7.8) 


(7.9) 


By  reason  of  its  dependence  on  (])  and  G‘^(r,F)  transforms  in  a  slightly 
more  complicated  fashion: 


G^(R(A(t)).r,F) 

^^(r,f(A(})).F) 


~  gjm((t) -(})•) 

m  =  -oo 

^  gjm((l)  -  0' ) 

m  =  -oo 


I 


ej'''*'l’'R(A«.G^(r,F)ej”^1> 


n  =  — oo 


I  e 
n  =  — oo 


inl^QS 

mn 


(r,F).R-l(A<l)) 
^j(nl^  -  m)A(t) 


(7.10) 

and 


G^(R(A())).r,R(A(|)).F)=  I 

m  =  -oo 

I  e-'"^^^'R(A(|)).G^^(r,F  ).I“^(A({))ej"^^^^ 
n  =  — oo 


(7.11) 
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7.3  The  dyadic  Green's  function  for  remote  source  and  field 
points 

Forr'— >00,  given  by  (3.37)  and  G^(r,F)  is  the  special  case  of  (7.3) 

with  the  currents  in  (7.4a)  corresponding  to  a  remote  source.  Substitution 

=s 

into  (7.6)  gives  the  elements  of  G^(f,F)  to  be 


8.  The  rotating  ensemble  in  the  presence  of  a  stationary 

ensemble 


8.1  The  ensembles 

We  now  consider  two  coaxial  ensembles  similar  to  that  studied  in  the 
foregoing  sections,  one  static  and  the  other  rotating  about  the  common  axis. 
This  is  illustrated  in  figure  8.1. 

The  stationary  ensemble,  say,  consists  of  P  bodies,  with  surfaces  Sj  , 
i=0,l,...I^-l,  the  ith  generated  from  the  0th  by  rotation  through  the  angle 

cD?=~i  .  (8.1) 

1 

The  rotating  ensemble,  say,  consists  of  bodies,  with  surfaces,  at  time  t  = 
0,  S^,  i=0,l,...I^-l,  the  ith  generated  from  the  0th  by  rotation  through  the 
angle 

.  (8.2) 

1  jR 

This  ensemble  rotates  about  the  z-axis  in  the  +({i  direction  with  angular 
velocity  Q. 

The  stationary  ensemble  is  the  environment  in  which  we  shall  study 
scattering  from  the  rotating  ensemble.  The  dyadic  Green's  function  for  this 

=s 

environment,  G  (r,F)  is  given  in  (7.5). 

Let  a  point  on  body  of  the  ensemble  at  reference  time  t  =  0  be  r^,  with 

azimuthal  angle  (j)..  At  time  t,  the  ensemble  has  rotated  through  angle  Qt, 

and  r.  and  (b.  have  become 
1 


r-  =  R(Qt).ri  =  R(<I>i  +Qt).ro 


(8.3a) 


({)|  =  (t)i  +  Qt  =  <1)0  +  Oi  +  Qt 


(8.3b) 


—t  Rt 

respectively,  and  r:  resides  on  a  surface  we  shall  call  .  We  extend  the 

concept  of  corresponding  points,  introduced  after  (2.3),  to  include  r|,  rj  and 

TQ  Whenever  similarly  sub/superscripted  position  vectors  or  their 
coordinates  appear  in  the  same  equation  they  are  to  be  understood  as 
corresponding. 
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8.2  Currents  on  the  rotating  ensemble 

The  incident  field  at  the  point  r| ,  due  to  a  dipole  of  moment  is  p  at  F  acting 
in  the  presence  of  the  stationary  ensemble,  is 

E'(r[)=G^(r|.F).p  (8.4) 

Let  the  current  density  distribution  at  rj  on  the  jth  body  be  jj(fj).  It  is 
convenient  to  define  a  related  function  J.(rQ,Qt),  through 

J](rJ)  =  R  +  Qt).  Jj(ro,Qt).  (8.5) 

Recognising  that  the  current  on  the  jth  body  at  time  t  is  the  same  as  that  on 
the  0th  at  the  time  when  it  occupies  the  same  position,  we  note 

Jj(ro.Qt)  =  Jo(ro,i^t  +  of ) 

J  (8.6) 

In  order  that  the  tangential  electric  field  on  the  rotating  ensemble  be  zero, 
Jj(rj)  must  satisfy  the  integral  equation  (Cf.(5.4)) 

1^-1= 

[  II  I  GS(f' fh.J!(fhdsV]  =-[G®(r!,F).pl 

f‘'onsR'j  =  0  JJ  J'  '  3 


J  J 


all  on  all  i 
1  1 


With  substitution  from  (7.5),  (7.7),  (7.8),  (8.3)  and  (8.5),  and  some 
rearrangement,  this  becomes 


1^-1  oo  oo  jm'((t). -({).■)  jn'I^(4).'+at) 


[  J  Z  I  I  e 

r^'cnS^  j  =  0m'  =  -<^n'=-oo 


1  J 


°°  jm((t)^ -(|)' )  jm(Qt  +  0?^)=q 

=  -[  I  e'  "0  .  GS(ro,F).p], 

m  =  —oo 

R 

all  Tq  on  Sq  ,  all  i 


where 
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,F).p 


^  jnI^(t)'7=rS  /- 
n  =  — oo 


(8.9) 


The  current  distributions  on  the  bodies  are  to  be  found  from  the  solutions  of 
the  above  integral  equation. 

Unlike  the  case  with  (5.8)  there  appear  to  be  no  simple  relations  between  the 
Jj(rQ',Qt)  from  one  body  to  another  or  over  time  or  between  directions  of 

incidence.  Nor  is  it  possible  to  reduce  the  surfaces  over  which  the  equation 

is  enforced  to  a  single  one.  To  solve  (8.8)  by  the  Moment  Method  it  will  be 
necessary  to  distribute  basis  functions  over  all  the  bodies,  and  enforce  the 
equation  at  a  number  of  points  equal  to  the  number  of  bases.  This  must  be 
done  for  a  series  of  time  values.  And  each  direction  of  incidence  requires  a 
separate  solution.  Thus  the  method  does  not  see  the  large  reduction  in  matrix 
dimension  enjoyed  by  the  earlier  equation  (5.8). 


The  solutions  Jj(rQ,^]!t)  are  periodic  with  fundamental  frequency  SI.  Their 

spectra  are  lines  separated  by  this  frequency.  It  is  sufficient  to  determine 
only  one  spectrum,  that  of  jQ(rQ,Qt)  say,  since  all  the  other  currents  are 

time-shifted  versions  of  this,  see  (8.6).  Let  the  number  of  lines  to  be 
determined  be  2N+1  (one  at  the  centre  and  N  on  each  side).  Then  N  is  the 
integer  which  most  closely  satisfies 

2N  +  1>”  (8.10) 


where  AQ.  is  the  frequency  bandwidth  of  jQ(rQ,Qt),  and  will  be  discussed 
below. 

Determination  of  the  spectrum  through  use  of  a  discrete  Fourier  transform 
requires  values  of  JoCrQ.i^t)  at  2N+1  times,  evenly  spaced  through  the 
2k  ^  2k 

rotation  period  —  at  angular  intervals  A(Qt)  = - .  It  is  noted  that  the 

n  2N-bl 

solution  of  (8.8)  for  a  given  ntjj  =nA(Qt),  yields  Jj(fo,ntn)  for  j  =  0...I^,  and 

-  n  j  p 

with  (8.6)  we  therefore  have  Jo(ro,(- - + -^)27t),  J  =  0...I  ;  i.e.  each 

2N-I-1 

solution  of  (8.8)  for  a  given  yields  jQ(rQ,Qt)  at  values  of  Qt.  For  2N+1 

spectral  components  we  therefore  need  to  solve  the  integral  equation 

XT  2N4-1 

Nsolve  =  times,  for  Qtn  =  nA(nt) ,  n  =  0,...,Nsolve  “  ^  This  corresponds 

to  solutions  at 


Qtrr  = 


n  27t 


^solve 


,  n  -  0,...,N5oive  ^ 


(8.11) 


i.e.  at  a  set  of  Ngolve  equi-spaced  intervals  in  the  inter-body  angle. 
With  (8.10),  Ngoive  is  the  integer  which  most  closely  satisfies 
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I^solve  - 


AQ 


(8.12) 


8.3 


The  current  spectra 


Assuming  that  (8.8)  has  been  solved  a  sufficient  number  of  times  and  a 
discrete  Fourier  transform  to  have  been  applied,  we  may  write  jQ(rQ,Qt)  as 
the  series  of  frequency  components 


J0(r0«^^t)=  X  (8.12) 


where  Jo^(ro) ,  ^  =  0,±1,±2,...,±N  is  the  discrete  Fourier  transform  of 
Jo(ro,^2tn)  as  determined  above,  and,  with  (8.5)  and  (8.6), 

-t  t  =  R  ^  -  j^(Ot+of) 

Jj(^])=R((D]^+^2t).  2  %(ro)e  J  (8.13) 

i=-oo 

It  was  seen  in  [1,2]  that,  for  a  rotating  ensemble  in  free  space,  the  current 
spectrum  extends  in  principle  over  an  infinite  band  of  frequencies,  but 
becomes  negligible  outside  a  finite  range.  The  extremes  of  the  significant 
spectrum  correspond  to  the  Doppler  shifts  associated  with  the  greatest  linear 
velocities  of  any  point  on  the  rotating  ensemble,  as  viewed  by  the  source. 

This  is  greatest  when  the  source  lies  in  the  plane  of  rotation,  in  which  case 
the  significant  frequencies  extend  kQRniax^  (where  Rmax  is  the  maximum 
radius  of  the  ensemble)  above  and  below  the  centre.  The  current 
components  with  the  extreme  frequencies,  due  to  the  rotation,  radiate  fields 

with  spectra  extending  over  similar  bands,  again  understandable  in  terms  of 
Doppler  shift.  The  backscattered  field  has  significant  frequency  components 
extending  2koRniax^  above  and  below  the  centre. 

In  the  composite  of  the  stationary  and  rotating  ensembles,  currents  in  each 
ensemble  radiate  fields  which  induce  currents  in  the  other  ensemble.  Thus  a 
frequency  component  in  one  ensemble,  because  of  its  motion  relative  to  the 

other,  induces  a  spectrum  of  currents  in  the  other,  and  vice  versa.  A  set  of 
current  frequencies  in  one  extending  over  a  band  ±koRmax^  about  the 
centre  induces  a  set  of  frequencies  in  the  other  over  a  band  ±2koRixiax^’ 
which  in  turn  induces  a  band  +3koRniax^  first,  and  so  on.  The 

resulting  bandwidth  is  infinite.  However  each  extension  of  the  band  occurs 

through  electromagnetic  interaction  between  the  currents  in  the  two 
ensembles.  If  the  mutual  coupling  is  weak,  as  it  may  be  if  the  ensembles  are 
separated  by  more  than  a  small  part  of  a  wavelength,  each  extension  is 
weaker  than  the  previous  by  a  constant  ratio.  The  current  spectrum  is 
therefore  expected  to  be  significant  over  the  first  few  ranges,  two  or  three 

times  IkoRinax^  perhaps,  and  to  decrease  exponentially  outside  this. 

An  estimate  of  AQ  is  therefore  2NintkoRmax^-  where  is  the  number  of 

interactions  between  the  ensembles  considered  significant  (two  or  three, 
say).  Then  with  (8.12),  and  noting  that  the  arc  distance  between  the  points 

27C 

on  the  bodies  at  maximum  radius  is  d  =  Rfnax~D">  we  have 
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^solve  -  ^intdX/2 


where  dj^/2  d  measured  in  half-wavelengths  of  the  illumination 
frequency . 


The  above  simple  relation  gives  the  number  of  times  (8.8)  must  be  solved. 
With  (8.11)  we  have  the  rotational  positions  at  which  solutions  are  required: 


Qtn  =  n 


27t; 


Nintdx/2 


,  n  =  0, 


Nintdx/2  ^ 


i.e.  at  each  of  NjfjtdX/2  equi-spaced  subdivisions  of  the  angle  between  the 
bodies. 

Experience  of  calculating  such  spectra  suggests  that  the  value  =  2 

provides  good  results  [7]. 

The  above  formulas  are  essentially  a  reflection  of  the  Nyquist  sampling 
criterion.  Seen  in  this  light  they  therefore  indicate  the  number  and 
separation  of  determinations  of  scattering  which  are  necessary  whether  by 
the  present  method  or  alternative  calculation  or  measurement. 
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8.4  The  scattered  field 

The  field  radiated  by  the  rotor  currents  is  the  sum  of  the  fields  radiated  by  all 
the  currents  on  all  the  bodies,  in  the  presence  of  the  stationary  ensemble 
Thus, 


E®(r,t)=  I  IJ  G^(f,r|').j|(rhds[’. 

i  =  0rt.,„sRt 


(8.14) 


With  (7.10),  and  noting  that  the  rotation  from  Tq'  to  f:'  involves  the 

D  =C _ f 

azimuthal  rotation  Atj)  =  Or  +  f2t ,  we  represent  G  (r,r/)  in  the  form 


G^(f.r|')=  S  I  e 


jm(({)-<t)o')gjnI^<j)0' 


m  =  — oon  =  — 


JCnl^  mXO? 


(8.15) 


With  (8.13)  and  (8.15),  (8.14)  becomes  after  rearrangement 


E^(r.t)=  six  gj(nlS-m+Oflgm(t) 

^=_ooni  =  -oon  =  — oo 

fo’onSo 

I-l  j(nI^-m+^)<Df 

X  e  J 


(8.16) 


From  the  definition  of  Ov^,  see  (8.2),  the  summation  on  i  is  zero  unless 

c  R  R 

nr  —  m  +  ^  is  an  integer  multiple,  h  say,  of  I  ,  in  which  case  it  equals  I  . 

R  S 

Therefore,  set  ^  =  hl  -nl  +m  and  sum  over  h.  (8.16)  now  takes  the  form 


E^(r,t)=  X  E^  „(f)e-'^“0'^‘^^ 

h  =  —oo  hi 


(8.17) 


where 


E*R(f)  =  lR  X  X  +  11  I 

hi  £  =  -oon  =  -oo  fo'onSQ'^  =  -~ 


^  S  R  .  (>'.'‘0')  Jo,^(«'0'))dso' 
nr-hl*^+^,n 


(8.18) 
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and  the  hitherto  suppressed  time  dependence 


eJ“0l 


has  been  reinserted. 
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8.5  The  spectrum  of  the  scattered  field 

(8.17)  represents  the  field  as  a  spectrum  of  components,  with  radian 
frequencies 


CO  „  =  (0„  +  hl^Q,  h  =  0,±1,±2 . ±°o  (8.19) 

hl^  ^ 


This  is  consistent  with  expectations.  In  time  —5 —  each  body  rotates  to  the 

position  previously  occupied  by  its  nearest  neighbour  and  the  ensemble  re¬ 
presents  the  same  scattering  geometry.  Thus  the  scattered  field  is  modulated 

periodically  with  the  fundamental  radian  frequency  I  Q  and  the  spectrum 
consists  of  lines  separated  by  this  frequency. 

An  insight  drawn  from  the  earlier  studies  [1,2]  is  that  the  values  of  indices 

C  D  =C  _ 

nl  -hi  +£  and  n  of  G  q  15  (r.fn')  are  limited  to  finite  ranges. 

nr-hl^-H^n 

However  it  has  been  seen  in  Section  8.3  above  that  the  current  index  £  may 
have  an  extended  range  due  to  multiple  interactions  between  the  stationary 

I) 

and  rotating  ensembles.  It  follows  that  hi  may  have  a  correspondingly 
extended  range. 

Thus  the  spectrum  of  the  rotating  ensemble  in  the  presence  of  the  stationary 
bodies  extends  beyond  the  limits  observed  in  the  case  of  the  ensemble 
rotating  in  free  space.  As  observed  for  the  induced  currents,  this  may  be 
attributed  to  the  electromagnetic  interaction  between  the  two  ensembles,  and 
will  be  greater  or  less  as  the  interaction  is  stronger  or  weaker.  It  may  be 
understood  as  a  Doppler  effect  due  to  the  changing  distances  between  the 
interacting  bodies. 

A  clear  example  of  the  extended  spectrum  is  seen  in  the  case  of  axial 
backscatter.  For  the  ensemble  rotating  in  free  space  there  are  no  frequency 
components  other  than  the  frequency  of  illumination.  In  the  presence  of 
the  stationary  ensemble,  the  spectrum  has  been  found  by  independent 
calculations  to  extend  over  a  range  of  components[7]. 
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Figure  8. 


Two  coaxial  ensembles,  the  upper  rotating,  the  lower  stationary. 
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9.  Conclusions 

The  study  of  electromagnetic  scatter  from  rotating  objects  with  angular 
periodicity  has  proceeded  through  three  stages:  (i)  modelling  of  the  objects 
as  wires  radiating  from  an  axis,  reported  in  [1];  (ii)  generalisation  of  the 
wires  to  include  bodies  of  arbitrary  shape,  reported  in  [2];  and  (iii)  the 
inclusion  of  an  ensemble  of  stationary  bodies  coaxial  with  the  rotating 
ensemble,  the  subject  of  the  present  report. 

A  procedure  has  been  developed  by  which  the  spectrum  of  the  scattered  field 
may  be  calculated.  As  has  been  the  case  previously,  it  is  possible  to  deduce 
some  of  the  characteristics  of  the  spectrum  from  the  mathematical  form  of 
the  equations,  without  proceeding  to  numerical  solution.  The  main 
conclusions  of  this  nature  are  that  the  spectral  lines  are  spaced  at  intervals 
of  the  "body  rate",  i.e.  the  rotation  frequency  multiplied  by  the  number  of 
bodies  in  the  rotating  ensemble,  and  that  the  spectrum  extends  beyond  the 
limits  of  the  Doppler  shift  associated  with  the  linear  velocity  of  the  bodies. 
These  conclusions  are  supported  by  independent  calculations  [7]. 

A  question  which  is  not  answered  here  concerns  the  possibility  of  specially 
strong  spectral  components  associated  with  simple  combinations  of  the 

R  S 

numbers  of  bodies  in  the  two  ensembles,  e.g.  (ml  +nl  )Q,  where  m  and  n 

S  R 

are  small  integers.  The  appearance  of  the  expressions  such  as  nl  —  hi  in 
(8.18)  is  intriguing,  but  no  conclusions  have  so  far  come  to  light.  If  such 
characteristics  exist,  it  may  be  necessary  to  seek  them  in  numerical 
solutions. 

It  would  be  possible  to  extend  these  studies  to  include  further  ensembles  of 
bodies,  stationary  and  rotating,  with  different  numbers  of  bodies.  Indeed,  in 
so  far  as  the  present  work  includes  the  extreme  case  where  there  is  only  one 
body  in  each  ensemble,  it  includes  multiple  layers  of  stationary  and 
synchronously  rotating  ensembles,  since  all  the  stationary  ensembles  may 
be  viewed  as  a  single  ensemble  of  one  body,  and  the  same  is  true  of  the 
synchronously  rotating  ensembles. 

However,  as  the  models  which  are  the  subjects  of  this  series  of  studies  have 
become  progressively  more  general  and  complicated,  the  opportunities  for 
exploiting  the  angular  periodicities  have  reduced.  The  analysis  has  become 
progressively  more  laborious  and  the  properties  of  the  spectrum  more 
difficult  to  state  categorically. 

For  these  reasons  it  seems  probable  that  further  progress  will  best  be  made 
through  numerical  and  experimental  studies.  The  proposed  continuation  of 
the  work  by  the  present  author  lies  in  the  first  of  these  directions. 
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Appendix  A 

The  following  formulas  may  be  found  in,  or  readily  derived  from  formulas 
found  in,  [6,  Ch  9]: 


Jpcos(t)^  Ijnj  (p)e-'"^ 


(Al) 


n  =  — oo 


cos 


n  =  — oo 

oo 


(A2) 


n  =  -oo 


(where  a  prime  denotes  differentiation  with  respect  to  the  argument) 
.  ,  jpcostt)  ” -n  n  jnd) 


n  =  — oo 

oo 


(A3) 


1  °° 

=  -r  i;j"(J„.i(P)  +  J„^,(p))e 


jn(|> 


n  =  — oo 


(A4) 


(A5) 
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